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SUMMARY 

An  analysis  was  made  of  flight  test  data  to  obtain  information 
representing  the  atmospheric  turbulence  inputs  to  the  test  aircraft.  The 
basis  of  the  method  was  the  use  of  the  linearized  equations  of  motion  of 
the  aircraft,  with  the  addition  of  certain  terms  representing  the  aerody- 
namic effects  of  the  turbulence  on  the  aircraft. 

These  equations  were  set  up  on  an  analog  computer  and  verified 
by  the  technique  of  matching  transient  response  flight  data.  Then,  through 
the  use  of  a  feedback  system,  the  desired  turbulence  quantities  were  ob- 
tained as  system  outputs.   Inputs  to  the  system  were  the  aircraft  respi-^nsa 
and  control  deflections,  obtained  from  the  flight  test  data. 

The  method  was  found  to  be  valid,  assuming  the  a'lrcraff  analogue 
was  correct.  The  resulting  data  appeared  suitable  for  further  analysis  to 
determine  Its  statistical  qualities. 
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I   INTRODUCTION 

In  accordance  with  a  proposal  submftted  to  the  Bureau  of  Naval 
Weapons,  the  James  Forrestal  Research  Center  at  Princeton  University  in 
late  1962  undertook  an  investigation  into  the  lateral  flying  qualities  of 
operational  Navy  aircraft  in  the  PA  configuration  approach  to  carrier 
landings. 

As  noted  in  Reference  I,  the  official   proposal,  it  is  typical 
that  during  this  type  of  approach  the  lateral  stability  and  control  charac- 
teristics of  some  high  performance  aircraft  deteriorate   with  increasi ng 
angle  of  attack  and  decreasing  airspeed.  This  is  common  in  many  current 
sweptwing  aircraft.  The  situation  is  aggravated  by  the  added  difficulty 
of  approach  and  landing  on  aircraft  carriers,  where  the  pilot  is  required 
to  maintain  a  f I i ghtpath  with  extremely  narrow  tolerances.  The  difficulty 
is  further  compounded  by  natural  atmospheric  disturbances,  and  those  caused 
by  the  airflow  about  the  carrier  hull  and  superstructure,  which  are  often 
encountered  in  the  final  few  seconds  prior  to  touch-down. 

It  was  therefore  proposed  that  a  flying  qualities  research  pro- 
gram be  conducted  at  Princeton,  utilizing  a  variable  stability  aircraft  to 
obtain  pilot  opinion  ratings.  The  approach  to  a  carrier  was  to  be  simula- 
ted as  closely  as  possible,  including  effects  of  the  forementioned  atmos- 
pheric disturbances.  These  turbulence  effects  were  to  take  the  form  of 
recorded  signals  applied  as  inputs  to  the  electronic  servos  driving  the 
control  surfaces  of  the  variable  stability  aircraft.   It  was  therefore  con- 
sidered desirable  to  represent  the  turbulence  as  accurately  as  possible 
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sfnce  this  factor  would  probably  have  an  Important  effect  on  pilot  rating 
in  this  particular  investigation. 

To  this  end,  a  suitably  instrumented  Navy  fighter  aircraft  from 
the  Naval  Air  Test  Center,  Patuxent  River,  Maryland,  made  a  series  of  ac- 
tual carrier  approaches  under  varying  conditions  of  wind  velocity  across 
the  deck  and  cross  wind  component.  Recordings  were  made  of  the  aircraft 
response  in  five  degrees  of  freedom,  together  with  recordings  of  the  var- 
ious control  surface  deflections.  These  data  were  analyzed  on  an  analog 
computer  setup  of  the  aircraft  equations  of  motion.  By  a  process  fo  be 
described,  it  was  possible  to  determine  the  atmospheric  disturbance  infoi — 
mation  in  a  form  suitable  for  application  to  the  variable  stability  aircraff. 

This  report  is  concerned  with  the  first  phase  of  this  research 
project:  the  acquisition,  reduction,  and  analysis  of  the  flight  test  data 
associated  with  the  carrier  approaches  made  by  the  Naval  Air   Test  Center 
aircraft. 
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II  EQUIPIChfT  AND  PROCEDURE 

I)  FLIGHT  TEST  EQUIPICNT 

The  data  acquisition  flights  were  flown  in  an  F4B  fighter  aii — 
craft  pi loted  by  Navy  test  pi  lots  from  the  Naval  Air  Test  Center.  The 
F46  (Phantom  II)  is  a  two-place,  swept-wing,  high-performance  interceptor 
manufactured  by  the  McDonnell  Aircraft  Corporation  of  St.  Louis.   It  is 
powered  by  two  General  Electric  J-79  turbojet  engines,  each  producing  a 
maximum  of  17,000  pounds  thrust  at  sea  level. 

The  aircraft  was  completely  instrumented  for  the  data  acqu/s/tlon 
by  the  Instrumentation  Branch  of  the  Flight  Test  Division,  Naval   Air   Test 
Center.  The  rear  cockpit  seat  was  removed  from  the  aircraft  to  provide 
space  for  the  instrumentation  padoge.  A  summary  of  the  flight  instrumen- 
tation fol lows: 

Power  Supply:  The  basic  aircraft  power  package  was  utilized  as 
the  instrumentation  power  source.  This  consisted  of  the  normal  28  volt 
DC  supply,  and  115  volt,  400  cycle  AC  source. 

Sensing  and  Transducing;  The  physical  quantities  measured  by  the 
Instrumentation  system  are  shown  in  Table  I.  Potentiometer  type  transducers 
were  chosen  in  order  to  provide  compatibility  with  the  telemetering  system. 
The  potentiometers  provided  signals  of  adequate  magnitude,  and  through 
their  excitation  source  also  provided  zero  and  full  scale  reference  vol- 
tage for  the  telemetering  unit. 

Table  I  also  lists  the  basic  instrument  used  in  the  measurement 
of  each  quantity.  The  rate  gyros  which  sensed  the  various  angular  velo- 
cities, and  the  potentiometers  which  measured  the  required  position  angles, 
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were  all  types  commonly  used  in  flight  test  Instrumentation,   It  should 
be  noted  that  the  actual  control  surface  deflections  were  measured  rather 
than  simple  stick  and  pedal  deflections,  due  to  the  non-linear  charactei — 
istics  of  the  stability  augmentation  system.  Use  of  the  actual  deflec- 
tions permitted  the  aircraft  equations  of  motion  to  be  written  in  simple 
form. 

Airspeed  and  altitude  were  measured  by  W,  0.  Leonard  force  balance 
type  transducers.  These  units  were  located  on  the  lower  shelf  of  the  rear 
cockpit  instrumentation  sled,  as  shown  In  Reference  2.  The  transducers, 
operating  on  the  force  balance  principle,  provided  linear  outputs  for  afi 
speed  and  altitude  without  the  use  of  cams,  function  gears,  shunts,  etc. 
The  force  sensors  were  Isolated  from  the  output  devices,  thus  minimizing 
response  to  transients.  As  a  result,  the  units  were  characterized  by  a 
pressure  accuracy  of  up  to  .\%,      In  the  case  of  the  altimeter  this  meant 
accuracy  to  within  three  Inches  for  a  two  hundred  foot  aircraft  altitude. 

Lateral  and  normal  acceleration  were  sensed  by  Conner  model  4310 
linear  accelerometers  Installed  with  a  range  of  +  .5  G  and  +  2.0  G,  respec- 
tively. These  units  exhibited  natural  frequencies  much  higher  than  the 
test  aircraft,  thus  precluding  resonance  problems.  Their  response  was  es- 
sentially flat  to  8  cps.   It  was  determined,  however,  during  Installation, 
that  some  H.  F.  system  noise  was  being  passed  by  the  accelerometer.   In 
order  to  eliminate  this  noise,  a  low  pass  filter  was  therefore  designed 
and  installed  between  the  accelerometer  output  and  the  telemetering  unit. 

The  filter,  designed  for  a  5  cps  corner  frequency,  eliminated 
the  high  frequency  noise  but  provided  essentially  flat  response  over  the 
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expected  aircraft  response  range,  which  was  estimated  to  be  no  higher  than 
I  cps,  v/ell  above  the  Dutch  roll  frequency, 

A  schematic  of  the  filter  circuit  is  shown  below,  together  with 
a  sketch  of  the  filter  frequency  response  characteristic: 
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Telemetering  and  Recording;   Telemetering  equipment  of  the  pulse 
width  type  was  installed  in  the  aircraft  as  part  of  the  instrumentation 
package.  This  equipment,  manufactured  by  the  Applied  Science  Corporation 
of  Princeton  (ASCOP),  was  the  manufacturer's  D  Series  Pulse  Width  Multi- 
coder-Telemeter,  capable  of  sampling  up  to  43  channels  at  a  sampling  rate 
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of  20  times  per  channel  per  second,  using  a  Model  DSG20  Commutator -Gear 
Box, 

Since  each  channel  was  sampled  20  times  per  second,  it  was  felt 
that  good  fidelity  would  be  obtained  with  input  frequencies  as  high  as 
the  2  to  5  cps  range.  Many  of  the  measured  quantities  were  sampled  on  each 
of  two  or  three  channels,  and  in  these  cases  it  was  felt  that  good  repro- 
duction was  assured  at  all  frequencies  encountered  in  the  investigation. 
Table  11  lists  the  measured  quantities  and  the  number  of  sampling  channels 
used  for  each.  For  those  quantities  sampled  in  more  than  one  channel,   the 
sampling  was,  as  far  as  possible,  equally  spaced  in  time. 

The  sampled  data  were  then  converted  into  PW  data  by  the  Model 
DKH-I  Keyer. 

The  pulse  width  information  was  then  recorded  directly  in  the 
aircraft  by  means  of  an  Ampex  Model  800  tape  transport  which  was  mounted, 
together  with  its  associated  power  supply,  in  the  rear  cockpit  installation. 

Provisions  were  also  made  for  simultaneous  voice  recording  by 
the  test  pi  lot. 

2)  GROUND  EQUIPMENT 

The  following  units  were  available  at  the  Forrestal  Research 
Center  for  analysis  of  the  flight  test  data: 

Magnetic  Tape  Recorder;  An  Ampex  Model  309C  dual  track,  three 
speed  tape  recorder  was  available  for  playback  of  the  recorded  flight  data, 
and  also  for  use  in  the  preliminary  analysis  studies  which  were  conducted. 

Telemeterino  Ground  Station:  The  telemetering  unit  used  for  de- 
coding the  taped  pulse  width  flight  test  data  was  the  ASCOP  M  Series  PW 
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Ground  Station,  This  unit  continuously  sampled  the  recorded  pulse  width 
information  and  converted  this  data  into  continuous  voltage  form.  The 
information  in  this  form  was  then  compatible  with  the  analog  computer  for 
analysis  purposes. 

GEDA  Analog  Computer:  The  basic  tool  in  the  analysis  process 
was  the  Model  L3  (GEDA)  linear  electronic  differential  analyzer,  manufac- 
tured by  the  Goodyear  Aircraft  Corporation.  This  computer  provided  twenty- 
four  automatically  stabilized  DC  computing  amplifiers,  each  with  an  open 
loop  gain  of  greater  than  5  x  lO''',  and  with  negligible  drift.  The  computer 
incorporated  an  automatic  error  indicator,  and  had  a  guaranteed  accuracy 
of  \%,     There  were  provisions  for  accurately  setting  computer  hoard   po- 
tentiometers using  a  special  null  indicator. 

Data  Recorder;  Visual  presentation  of  data  was  made  possible 
by  a  Sanborn  Model  154-IOOB  four  channel  recorder.  This  instrument  had 
a  very  high  natural  frequency  (42  cps),  and  low  frequency  response  flat 
to  zero  cps. 

3)   INSTRUMENT  CALIBRATION 

Calibration  of  the  flight  test  instrumentation  was  conducted  by 
the  Instrumentation  Branch  of  the  Flight  Test  Division,  Naval  Air  Test 
Center.  The  calibration  charts  are  shown  in  Reference  2.  All  gradients 
were  determined  to  be  essentially  linear  within  the  range  of  measured  data 
encountered  during  the  test  flights. 
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4)  PROCEDURE 

Data  Acqui  sition ! 

In  December,  1962,  the  Instrumented  F4B  aircraft,  piloted  by  a 
Navy  test  pilot,  made  a  series  of  thirty-three  landing  approaches  to  the 
aircraft  carrier  USS  Enterprise.  These  approaches  and  landings  were  made 
over  a  span  of  five  separate  flights.  The  heading  and  speed  of  the  cai — 
rier  were  varied  during  the  course  of  the  series  to  provide  a  number  of 
different  conditions  of  wind  velocity  across  the  deck,  and  cross  wind  com- 
ponent. The  wind  conditions  for  each  approach  are  shown  in  Table  111. 

It  was  expected  that  the  most  severe  turbulence  conditions  would 
be  experienced  in  the  presence  of  a  "right  to  left"  crosswind.  This  was 
due  to  the  extremely  turbulent  flow  about  the  carrier  "island"  being  direc- 
ted Into  the  path  of  an  approaching  aircraft,  usually  at  the  critical  point 
occurring  a  few  seconds  prior  to  touchdown.   It  was  felt  that  by  making 
the  approaches  under  varying  wind  conditions,  it  would  be  possible  to  ob- 
tain turbulence  data  which  would  be  representative  of  the  full  range  of 
severity  normally  encountered.  As  noted  in  Reference  3,  the  turbulence 
"burble"  associated  with  the  Enterprise  was  probably  the  extreme,  due  to 
the  very  box  I  Ike  design  of  the  ship's  "Island". 

The  Naval  Air  Test  Center  provided  facilities  aboard  the  Entei — 
prise  for  tape  playback  after  each  test  flight.  On  the  spot  assessment  of 
data  quality  was  thus  possible,  and  minor   adjustments  v/ere  made,    as  neces- 
sary, in  the  instrumentation. 
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Aircraft  Calibration: 

Since  the  basis  of  a  successful  analysis  of  the  flight  test  data 
was  to  be  the  setup  of  the  aircraft  equations  of  motion  on  the  analog  com- 
puter, it  was  necessary,  as  the  next  step,  to  obtain  a  valid  set  of  equa- 
tions representing  the  F4B  aircraft  in  the  PA  configuration.  The  basic 
method  used  was  the  matching,  on  the  computer,  of  transient  response  flight 
test  data.  The  method  will  be  discussed  in  some  detail  in  a  subsequent 
section. 

To  this  end,  two  additional  flights  were  made  of  the  test  air- 
craft for  the  purpose  of  obtaining  transient  response  data.  The  aircraft 
was  flown  in  the  PA  configuration,  and  pilot  control  inputs  of  the  pulse, 
step,  and  doublet  types  were  made.  The  aircraft  response  was  recorded 
as  before. 

In  order  to  prevent  stability  augmentation  equipment  inputs  to 
the  system,  these  units  were  made  inoperative  during  the  transient  res- 
ponse data  runs. 
Data  Analvsi  si 

Having  thus  acquired  the  necessary  flight  test  data,  the  actual 
analysis  was  then  conducted  at  the  Forrestal  Research  Center,  This  analysis 
is  discussed  in  detail  in  the  following  section  of  this  report. 
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1 1 1  DATA  ANALYS IS:  DESCR  f  FT  I  ON  AND  D I SCUSS I  ON 

The  first  step  in  the  data  analysis  was  the  determination  of  a 
valid  set  of  differential  equations  representing  the  F4B  aircraft  In  the 
PA  configuration.   It  was  assumed  that  the  motions  of  the  aircraft  while 
approaching  a  carrier  in  this  configuration  could  be  described  by  the  usual 
linear  equations.  This  assumption  seemed  valid  since  the  final  part  of 
a  typical  carrier  approach,  where  the  turbulence  effects  are  normally  most 
severe,  is  basically  a  constant  power,  constant  airspeed,  shallow  descent 
to  a  touchdown;  the  aircraft  motions  are  approximately  perturbations  from 
the  steady  state  glide  path. 

The  resulting  equations  are  shown  below,  written  in  dimensional 
form  for  compatibility  with  the  flight  data  and  analog  computer.  For  con- 
venience, the  terms  of  each  equation  have  been  divided  by  the  coefficient 
of  the  highest  appropriate  derivative.  The  assumption  was  also  made  that 
the  aircraft  forward  velocity  remained  essentially  constant  during  the 
runs,  thus  eliminating  the  drag  equation,  and  also  eliminating  longitudinal 
velocity  terms  from  the  remaining  equations, 

Lonqitudi  nal ! 

(1)  lift  equation:   a  -  3  +  [  ^h^  ]  a  +  [  ^lis  ]  65  =  0 
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(5)  roll    equation:    ^-C  %_.  ]c^  -   [  f^  ]t  -  [  ^f;^  ]3   -  [  ^^-^^a   je^ 

2JxT  2Jj^t  J^t^  Jj^  t^ 


(6)     van  equation:     if  -  [  ^JW_  ]*  -  [  12p_  3$  -  C  !!_§.  ]P  "  [       "^a   ]6 

2J,T  2J,T  J,T^  J,   -^ 


^7^ 


These  uncoupled  sets  of  equations  were  then  separately  mechan- 
ized on  the  GEDA  analog  computer.  Schematic  circuit  diagrams  of  the  two 
setups  are  shown  In  Figures  I  and  2. 

All  factors  in  these  equations  were  known  accurately,  except 
for  the  various  stability  derivatives.  Approximate  values  for  the  deriva- 
tives, for  the  PA  configuration,  were  obtained  from  a  report  of  the 
McDonnell  Aircraft  Corporation.  These  values  determined  the  initial  po- 
tentiometer settings  for  the  computer  setups.  Shown  below  Is  a  list  of 
the  various  quantities  utilized  In  tiie  initial  formulation  of  the  aircraft 
equations  for  the  PA  configuration: 
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Vo 

= 

134  knots  =  224  feet  per  second 

w 

= 

34000  lbs. 

do 

= 

16  degrees 

Clo 

= 

I.I 

P 

= 

.002378  slug/ ft® 

b 

= 

39  feet 

S 

= 

530  feet^ 

c 

= 

16  feet 

Ix 

= 

23880  slug  -  ft^ 

'y 

= 

121007  slug  -  ft^ 

Iz 

= 

131 177  slug  -  ft^ 

C,   =  3.05/radlan 
a 

Cm   =  --'9 
"a 


C.    =  +.328 
^s 


C.   =  -.240 

C^^   =  +.055 

C.   -  -.010 

Op 

Cff        =  -.285 
^P 
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C.  =  +.225 

Cnp  =  +.205 

Cnsg  =  -.006 

Cnp  =  +.036 

Cn^  =  -.330 

S  =  -.625 

Cv  =  +.950 


Cw   =  +.  780 


The  equations  of  motion  based  on  these  values  are   as  foJJows: 
Longitudinal' 

1.  a  -  9  +  .406  a  +  .0436  65  =  0 

2.  e  +  .  162  9  +  .067  a  +  .357  a  +  .792  63  =  0 

3.  n  +  2.81  a  -  .311  6s  =  0 

Latera I ; 

4.  -  p  -  .0832  3  -  .990  Jf  +  .0109  cp  +  .  139  cp  +  .0166  6r  =  0 

5.  -9  -  1.25  (p  +  .986  i  -  12.2  p  -  .510  6p  +  2.80  6^=0 

6.  -  J  -  .264  }  +  .0288  cp  +  1.89  p  -  .693  6^  -  .0556  63  =  0 

Improving  the  validity  of  these  equations  as  representations  of 
the  aircraft  thus  requirod  a  more  acriir<3+e   determination  of   the  stability 
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derivatfves.  This  was  accomplished  by  a  transient  response  curve  fitting 
technique,  such  as  outlined  In  Reference  A,   and  described  in  detail  In 
Reference  5,  The  mathematical  basis  for  the  validity  of  the  method  is 
found  in  Reference  6, 

For  this  method  a  number  of  the  transient  response  runs  were 
selected  by  reference  to  Sanborn  recordings  of  the  aircraft  response.  The 
selection  was  made  on  the  basis  of  freedom  from  telemetering  dropouts,  and 
relative  freedom  from  apparent  extraneous  inputs,  such  as  turbulence  effects, 
These  response  curves  were  of  the  rudder  pulse,  rudder  doublet,  aileron 
step,  aileron  doublet,  stabi  lator  pulse,  and  stabilator  doublet  types. 

Transparent  overlays  were  made  from  the  Sanborn  records,  incoi — 
porating  the  scale  changes  necessary  to  provide  agreement  with  the  scaling 
of  the  computer  outputs.  It  was  then  possible,  by  means  of  the  ground 
station  equipment,  to  apply  voltages  to  the  aircraft  analogue  represent- 
ing the  actual  aircraft  control  motions,  as  forcing  functions  to  either  the 
longitudinal  or  lateral  sets,  as  appropriate.  Sanborn  recordings  of  the 
analog  responses  were  made,  all  to  the  same  scale  as  the  transparent  over- 
lays. The  two  sets  of  data  were  compared,  and  by  an  iteration  procedure, 
the  potentiometer  settings  representing  the  various  stability  derivatives 
were  changed  until  close  agreement  was  attained  between  flight  response 
and  analog  response. 

In  the  lateral  case  the  matched  quantities  were  roll  rate  and 
yaw  rate.  It  was  found  that  a  close  match  to  the  flight  data  could  be 
fairly  rapidly  obtained  by  corrections  to  those  derivatives  which  had  major 
effects  on  the  important  lateral  response  parameters,  i.e.,  Dutch  roll 
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perlod  and  dampfng,  initial  roll  and  yaw  magnitudes,  spiral  stability. 
Improvements  to  the  Initial  match  were  obtained  by  a  systematic  iteration 
procedure,  varying  the  potentiometer  settings  In  turn  until  the  desired 
match  resulted. 

In  the  Initial  study  of  the  lateral  responses  it  was  immediately 
apparent  that  the  flight  data  was  not  Ideal,  During  even  the  best  runs, 
the  aircraft  obviously  had  not  been  In  steady  flight  at  the  Instant  of  the 
control  Inputs;  there  were  initial  conditions  which  had  to  be  considered 
in  the  matching  process.  By  means  of  the  analog  initial  condition  cir- 
cuitry It  was  possible  to  apply  to  the  computer,  as  Initial  conditions, 
the  actual  flight  data  response  in  roll  rate  and  yaw  rate  existing  at  the 
instant  of  control  deflection.  This  was  considered  to  be  an  improvement 
but  it  was  noted  that  since  the  lateral  equations  represented  a  fourth 
order  system,  two  additional  initial  conditions,  3  and  cp,  for  example, 
were  necessary  for  accurate  solution.  These  were  not  available  and  some 
error  was  thus  to  be  expected.   In  the  Iteration  procedure  It  was  found 
that  the  best  matches  were  obtained  with  flight  data  where  the  aircraft 
was  very  close  to  steady  state  flight  conditions  at  the  time  of  control 
deflection.  This  was  to  be  reasonably  expected  since  In  these  cases  all 
Initial  conditions  were  approximately  zero. 

It  was  also  apparent  from  observation  of  the  flight  data  that 
some  atmospheric  turbulence  existed  during  the  tests.  This,  of  course, 
could  not  be  accounted  for  in  the  matching  process,  and  was  thus  a  source 
of  error.  For  this  reason  the  smoothest  runs  were  selected  for  the  match- 
ing process. 
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It  was  noted  also  that  during  the  various  transient  response  runs 
there  were  inadvertent  Inputs  of  those  control  surfaces  which  theoretically 
were  being  held  in  the  neutral  position.  These  inputs  were  also  played  into 
the  analog. 

It  can  be  seen  from  Table  I  that  only  the  left  hand  aileron  and 
spoiler  were  instrumented  for  deflection  measurement.  Since  the  equations 
of  motion  as  formulated  on  the  computer  were  written  in  terms  of  aileron 
inputs,  which  are  the  sum  of  right  and  left  aileron  deflection,  it  was 
necessary  to  derive  some  function  of  left  spoiler  deflection  which  might 
be  utilized  as  equivalent  right  aileron  deflection.  To  this  end,  a  cali- 
bration run  was  made  to  determine  the  relationship  between  aileron  and  spoiler 
deflection.  The  results  were  of  the  form  shown  in  Figure  3,  The  relation- 
ship between  the  two  was  shown  to  be  approximately  linear.  Therefore,  it 
was  possible  to  express  right  aileron  as  a  linear  function  of  the  recorded 
spoiler  deflection.  The  aileron  signal  fed  to  the  analog  was: 

a     a^    a^      aL      sp,_ 

It  can  be  seen  from  Figure  3  that  each  aileron  exhibited  a  small 
upward  deflection  range  of  approximately  three  degrees.  This  property  was 
ignored  in  determining  the  proportionality  constant  between  aileron  and 
spoiler.   It  was  felt  that  the  upward  deflections  were  primarily  for  the 
purpose  of  assisting  the  effect  of  the  spoiler  deflections  on  the  same 
wing.   In  the  slightly  open  position  the  spoilers  were  assumed  to  be  rela- 
tively ineffective  in  disrupting  the  airflow,  and  the  small  upward  aileron 
deflection  made  up  for  the  deficiency,  thus  maintaining  system  linearity. 
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By  means  of  the  analog  If  miter  circuits,  the  voltages  representing  these 
upward  aileron  deflections  were  eliminated  from  the  computer  inputs. 

It  should  be  noted  that  in  spite  of  the  turbulence  and  initial 
condition  problems,  it  was  possible  to  obtain  good  matches  between  flight 
data  and  analog  output.  The  accurate  determination  of  all  the  various 
stability  derivatives  would  probably  have  been  difficult  because  of  these 
problems,  but  the  primary  interest  here  was  the  matching  of  aircraft  re- 
sponse. This  could  be  done  by  matching  the    imporfanf    lateral  parameters 
mentioned  previously.  Since  these  parameters  are  primarily  control  led   by 
relatively  few  of  the  derivatives,  it  was  possible  to  obtain  a  good   match 
even  though  some  of  the  derivatives  might  have  been  in  error. 

It  should  also  be  noted  that  since  the  instrumentation  used  in 
obtaining  data  for  this  matching  was  the  same  as  subsequently  used  during 
the  actual  carrier  approaches,  then  to  the  extent  that  any  instrumentation 
errors  would  have  caused  false  matching,  the  effect  would  tend  to  be  can- 
celled during  the  actual  approach  data  acquisition. 

Shown  below  are  the  circuits  utilized  to  provide  inputs  to  the 
analog  from  the  recorded  data.  By  means  of  appropriate  potentiometer  set- 
tings the  quantities  were  all  scaled  correctly  for  compatibility  with  the 
analog  setup.  Biasing  of  all  quantities  to  the  correct  zero  position  was 
done  by  adjusting  the  biasing  knobs  of  the  ground  station  while  referring 
to  voltmeter  readings  of  the  various  tape  outputs. 

Shown  in  Figures  4,  5,  and  6  are  some  of  the  matched  transient 
response  curves.   It  is  to  be  noted  that  the  aileron  step  Inputs  were  used 
primarily  to  match  initial  roll  rate  amplitude.   It  was  felt  that  the  aircraft 
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was  out  of  the  linear  range  shortly  after  application  of  an  aileron  step 
and  this  precluded  any  additional  useful  matching. 
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The  longitudinal  responses  were  matched  by  comparing  pitch  rate 
and  normal  acceleration.  The  method  utilized  and  the  problems  encountered 
were  similar  to  those  discussed  in  reference  to  the  lateral  equations. 
The  phugoid  mode  was  not  visibly  excited  in  any  of  the  transient  runs  and, 
therefore,  the  matching  was  confined  solely  to  the  short  period  oscil- 
lations. The  matching  iteration  procedure  converged  much  more  rapidly  than 
in  the  lateral  case.  This  was  to  be  expected  since  in  this  longitudinal 
approximation  there  were  only  two  degrees  of  freedom  and,  consequently, 
fewer  controlling  derivatives  to  be  adjusted. 
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Shown  below  are  the  Input  circuits  to  the  longitudinal  analog: 


ffu/ht 
Ddta 


3C^  Una 


--^^— Otw^ 


pots 


^  ^^ 


e 


CO {^ \AA.'^ 


AAA/^ — I 


to  IC  'I 


In  Figures  7,  8,  and  9  are  shown  some  of  the  matched  transient 
response  curves. 

It  should  be  noted  that  the  Initial  matching,  both  longitudinally 
and  laterally,  showed  a  number  of  discrepancies  between  the  actual  aircraft 
response  and  the  analog  response  which  was  based  on  the  stability  deriva- 
tives obtained  from  the  McDonnell  report.   In  the  lateral  case  the  analog 
response  showed  a  predicted  Dutch  roll  period  of  approximately  4.5  seconds, 
with  neutral  damping.  As  can  be  seen  in  Figures  4  and  5,  the  aircraft 
actually  exhibited  a  well  damped  Dutch  roll  of  about  a  3  second  period. 
Longitudinally,  the  major  discrepancy  was  In  the  apparent  degree  of  air- 
craft static  stability,  which  was  considerably  greater  than  the  predicted 
value. 

These  discrepancies  are  illustrated  more  clearly  by  reference  to 
the  final  equations  of  motion  which  resulted  from  the  matching  process.   If 
these  equations,  which  are  shown  below,  are  compared  term  by  term  with  the 
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original  set,  the  changes  which  had  to  be  made  during  the  matching  are 
apparent. 

Lonqi tudi  nal  i 

(8)  a  -  9  +  .472  a  +  .0436  65  =  0 

(9)  e  +  .838  e  +  .067  a  +  9.00  a  +  4.50  65  =  0 

(10)  n  +  3.25  a  -  .311  6s  =  0 

Latera I ; 

(11)  -  3  -  .0956  p  -  .990  J  +  .0110  cp  +  .140  cp  +  .0209  6r  =  0 

(12)  -  9  -  2.48  cp  -  .984  jr  -  15.8  p  +  .497  dp  +  5.48  63  =  0 

(13)  -  \|r  -  .569  i  +  .0288  cp  +  3.38  p  -  .674  6r  -  .061  63  =  0 

It  can  be  seen  that  in  the  lateral  set  the  directional  stability 
and  yaw  damping,  which  are  the  primary  parameters  affecting  Dutch  roll 
period  and  damping,  respectively,  were  each  approximately  doubled  in  the 
matching  process.  Longitudinally,  it  can  be  seen  that  the  static  stability 
was  also  increased  considerably  during  the  matching.  Although  the  high 
value  of  Cfiu,/  which  finally  produced  a  good  match,  was  probably  unrealistic 
it  at  least  indicated  that  the  predicted  value  was  too  low.   It  is  noted 
again  that  accurate  determination  of  each  stability  derivative  was  by  no 
means  assured.  However,  the  aircraft  response  curves  were  fairly  accurately 
matched,  and  this  was  the  important  point  for  purposes  of  this  investigation. 

With  the  equations  of  motion  of  the  aircraft  thus  formulated, 
the  next  step  was  the  actual  analysis  of  the  flight  data  to  obtain  the 
turbulence  information.  Preliminary  studies  were  made  to  determine  the 
feasibility  of  the  proposed  method  of  solution.  For  these  studies  the 
lateral  equations  of  the  Navion  aircraft,  as  obtained  from  Reference  5, 
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were  set  up  on  the  analog  computer.  The  basis  of  the  method  was  the  well' 
known  feedback  principle: 


R 


^^ 


B 


C 


From  the  above  diagram,  C/R  =  A/(I  +  AB).   If  A,  the  forward 
loop  gain,  is  very  much  greater  than  unity,  then. 


C 
R 


B 


If  B  were  an  aircraft  transfer  function,  e.g.,  9/63^  then, 

C  «  _J pa    ia 

R     c?)/6a     cf) 

It  seemed  possible,  then,  that  if  known  aircraft  response  in  roll  rate,  for 
example,  were  fed  into  such  a  circuit  at  "R",  where  "B"  was  the  aircraft 
transfer  function  as  set  up  on  the  analog  computer,  then  the  circuit  out- 
put at  "C"  would  be  a  voltage  representing  the  aileron  deflection  which 
had  caused  the  roll  response. 

The  validity  of  this  theory  was  tested  by  recording  on  tape  the 
analog  response  i  n  cp,  ^,   and  p,  to  an  arbitrary  sequence  of  voltages  repre- 
senting aileron  Inputs  to  the  system.  This  was  accomplished  by  connect- 
ing a  potentiometer  to  the  appropriate  points  of  the  analog  circuit  and, 
in  effect,  "flying"  the  aircraft  in  the  lateral  modes.  The  previously 
discussed  feedback  circuit  was  then  connected  as  shown  below: 
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The  taped  response  in  roll  rate  was  then  played  into  the  circuit  at  "R" 
and  the  output  "C"  was  compared  with  the  original  aileron  signals.  The 
match  was  excellent  with  forward  loop  gains  as  low  as  10.  Below  this  value 
the  reproduction  decreased  in  accuracy  since  (C/R)  ?»  (I/B)  was  no  longer 
a  good  approximation. 

The  next  logical  step  was  to  study  the  effects  of  two  control 
inputs,  both  aileron  and  rudder.  Again  the  analog  was  "flown",  this  time 
with  aileron  and  rudder.  The  same  aircraft  response  quantities  were  again 
taped.  The  analysis  circuit  was  similar  to  the  previous  case,  with  the 
addition  of  another  loop.  The  circuit  Is  shown  below  in  block  diagram 
form.  The  taped  responses  in  9  and  ^  were  played  into  the  circuit  at  R 
and  R'  respectively,  and  the  outputs  at  C  and  C  were  compared  with  the 
original  69  and  6p.  Again  the  n)ei'ch   was.  evcellent,  down  to  forward  loop 
gains  of  10, 
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Flgure  10  shows  a  more  detailed  schematic  of  this  last  analysis 
circuit.   It  was  important  to  provide  negative  feedback  in  the  loops  in 
order  to  maintain  system  stability. 

The  results  of  this  preliminary  study  pointed  up  the  uniqueness 
of  the  solutions  to  the  aircraft  equations  of  motion.  A  given  combination 
of  aileron  and  rudder  inputs  results  in  a  unique  combination  of  roll,  yaw, 
and  sideslip  response. 

Originally,  it  had  been  proposed  that  this  exact  analysis  method 
be  employed  in  the  study  of  the  flight  test  data.   In  the  lateral  case, 
for  example,  aircraft  response  in  roll  rate  and  yaw  rate  were  to  be  applied 
to  a  similar  analysis  circuit,  where  the  analog  would  be  of  the  F4B.  From 
the  resulting  outputs  of  6^-  and  63  then  would  have  been  subtracted  the 
recorded  flight  test  control  deflections.  The  difference  would  be  an 
aileron  and  rudder  signal  i-epresent»ng  the  iurbulence   input  to  the  FAB, 


-24- 


Th  is  Information,  suitably  modified  to  account  for  the  difference  in  con- 
trol power  between  F4B  and  Navlon,  might  then  be  applied  to  the  variable 
stability  aircraft  as  a  simulation  of  the  turbulence  Input, 

The  longitudinal  analysis  would  be  similar.  Aircraft  response 
in  pitch  rate  would  be  applied  to  the  appropriate  circuit  and  the  output 
6g  would  yield  an  equivalent  65  representing  the  turbulence  input. 

This  method,  however,  obviously  neglects  the  linear  perturbations 
of  aircraft  motion  due  to  turbulence,  those  motions  along  the  aircraft 
axes.   It  assumes  that  the  atmospheric  turbulence  effects  may  be  des- 
cribed accurately  by  the  resulting  moments  about  the  aircraft  axes.  For 
purposes  of  this  research  project  reasonably  valid  simulation  may  be  at- 
tained based  on  this  concept. 

However,  In  the  interest  of  more  accurate  gust  simulation,  it  was 
decided  to  conduct  some  initial  research  into  the  possibility  of  obtaining 
turbulence  data  which  might  be  expressed  in  five  degrees  of  freedom,  thus 
including  any  lateral  and  normal  acceleration  effects.   In  addition,  the 
possibility  existed  that  some  qualitative  information  concerning  the  char- 
acter of  the  turbulence  itself  might  result  from  such  a  study. 

It  was  known  that  it  is  possible  to  visualize  the  atmosphere 
through  which  the  aircraft  moves  as  being  in  motion  at  some  steady  mean 
value,  with  random  turbulent  fluctuations  superimposed.  These  random  velo- 
cities may  be  considered  as  inputs  to  the  aircraft  linear  system,  which 
may  respond  in  six  degrees  of  freedom.  The  velocity  Inputs  are  assumed 
small,  permitting  linearization  of  the  resulting  aerodynamic  forces. 
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As  outlined  in  Reference  7,  it  is  convenient  to  introduce  these 
inputs  into  the  equations  of  motion  in  the  form  of  velocity  components 
Ug,  V  ,  and  Wg  (subscript  "g"  standing  for  gust).  If  it  is  assumed  that 
the  gust  wavelengths  are  all  large  relative  to  the  aircraft  dimensions, 
then  the  gradient,  or  velocity  variation,  of  these  gusts  across  the  length 
and  span  of  the  aircraft  may  be  considered  essentially  linear.  It  can 
be  shown  that  under  this  assumption  the  effects  of  the  gust  components  and 
their  gradients  are  aerodynamical ly  equivalent  to  perturbations  in  angle 
of  attack,  sideslip,  roll  rate,  yaw  rate,  and  pitch  rate.  These  are  noted 
as  ag,  3g,  9g,  ig,  Qg. 

The  aircraft  equations  of  motion  incorporating  the  modified  var- 
iables thus  become: 

Longitudinal ; 

•    P  C 

( 1 4)  a  -  9  +  Jia  (a  +  a„)   +  ""^s  6  s  =0 

2  T      ^    2t 

(15)  9  -  2llde  (9  +  9g)  -  ^  a  -  £lla  (a  +  aq)  -  ^  65=0 

Ht  Ht      Ht^      ^    Ht^ 

Latera 1 ; 

•   C  r  ^     r 

(16)  3  -  Ji.  (p  +  3q)  +  t  -  -In  (i  +  in^  +  -A  9  -  -Ie.  ^9  +  ^g^ 

2t      ^       4m,       ^2t     4|j, 

-  ^ySa  6,  -  ^yV  6,  =  0 
2t       2t 
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(17)  9  --^  (9  +cpq)  -^  (i  +L)  -!fi|(3  +3.)  -iffkSa 


JvT 


5-  -^ 


(18)  i-£[iL.(i'  + V  -S_(cp+cp,)  -!S  o  +3.)  -!^6, 

2J^T        ^     2J^T       ^     J^T^         ^     ^17^ 


J,T- 


6p  =  0 


With  the  gust  terms  rearranged  to  point   out  their  role   as  forcing  functions, 
the  equations  are  written  as: 
Lonqitudi  nal : 

(19)  d  -  0  +2^a  +^li^  6s  +^an  =0 

2  T     2t      2  t 

(20)  e  -  ^  §  -  ^  ji  -  £na  a  -  £2£^  6s  -  £lld£  §g  "  ^l^  ^g     =0 

Ht  Ht  Ht^  Ht^  Ht  Ht^ 

Lateral : 

(21)  3-!!i3+t-2l!:}+^cp-^cp-^6a-£y«r6p-£m3g 

2  T  4  ij.  2t  4  p.  2t  2t  2  t 


^Vr   ;         Cy 


4  |i  4  |j, 


(22)     $-%.<p-2f^t-!ff£3-!ff^63-''V6,-^*g 

2JxT  2JxT  JxT^  J^T^  JxT^  ZJ^T 

2JvT      3         J?^ 
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(23)     'i  -Sl^'^  '^^>  -^^  -  f2h.  6a  -  ''^"^r  6.  -  ^nr     I 

2J2T  2JzT  JzT^  JzT^  JzT^  2J2T       ^ 

It  I s  to  be  noted  that  these  equations  are  not  valid  for  short 
wavelength  gust  components. 

However,  It  can  be  shown  that  the  effect  on  aircraft  response  of 
the  higher  frequency  gust  components  is  negligible.  For  example,  the  ef- 
fective gust  frequency,  Og^  may  be  expressed  by: 

^g   =Vo. 
\ 

where  \   is  the  gust  wavelength.  For  a  short  X,  say  equivalent  to  the  air- 
craft wingspan,  Wg  I s  on  the  order  of  5  cps,  for  the  F4B  in  the  PA  configu- 
ration. This  is  very  much  higher  than  the  aircraft  natural  oscillatory 
frequencies  and  consequently  the  aircraft  response  is  negligible. 

In  addition,  for  gust  wavelengths  of  this  order,  and  shorter, 
the  aircraft  begins  to  "average  out"  the  gust  velocities,  with  the  result- 
ing net  gust  acting  on  the  airframe  being  near  zero.  Thus,  for  these  higher 
frequency  components  there  is  a  compounded  attenuation  -  that  due  to  the 
aircraft  dynamic  characteristics,  and  that  due  to  the  cancelling  effect 
just  described. 


r/%™ 
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|t  should  be  noted  that  if  the  F4B  were  to  respond  dynamically 
to  the  higher  frequencies,  this  last  effect,  that  of  cancelling,  would 
cause  some  error  if  the  acquired  turbulence  information  were  applied  to  the 
study  of  an  aircraft  differing  greatly  in  size  from  the  F4B.   It  is  quite 
possible  that  the  aircraft  in  question  might  respond  to  those  gusts  which 
were  cancelled  by  the  F4B  and  were  thus  missing  from  the  acquired  data. 
However,  all  aircraft  which  will  be  simulated  by  the  variable  stability 
Navion,  for  this  investigation,  are  similar  in  size  to  the  F4B,  and  there- 
fore could  be  expected  to  "average"  the  gusts  in  the  same  manner. 

With  these  modified  equations  of  motion  In  mind,  it  was  desired 
to  obtain,  as  a  solution,  the  various  gust  quantities.  The  method,  althouyh 
based  on  the  previously  discussed  feedback  principle,  can  probably  be  best 
described  by  a  single  degree  of  freedom  example.  For  this  purpose  it  is 
assumed  that  the  aircraft  can  respond  only  In  roll,  and  that  the  only  gust 
quantity  present  is  9  ,  The  equation  of  this  system  is  then: 

(24)  ^'  -  %_  cp  -  Ifiia  6a  -  ^^^^6r  6r  -  2iP_  cp_  =  0 

This  equation  is  set  up  on  the  analog  computer  and  a  feedback 
circuit  is  connected  as  shown  on  page  29. 

The  output  of  the  forward  loop  amplifier  Is  defined  as  -  cpg, 
and  this  quantity,  modified  by  the  appropriate  proportionality  constant, 
is  used  as  an  analog  forcing  function.  The  recorded  values  of  -  6p  and 
-  6a  are  also  used  to  force  the  computer  setup.  The  following  equation, 
which  expresses  the  aircraft  response.  Is  satisfied  at  point  A: 

^recorded  ~  "^control  response  ■*"  ^gust  response 
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0r<i(.orclecl   /  A 


< 


'i^conitrol  rasponst 


Assuming  that  the  circuit  is  dynamically  stable,  then  \he   -  (^g 
quantity  actually  appears  at  the  amplifier  output.  This  circuitry  is  simi- 
lar to  those  previously  discussed. 

The  more  complicated  circuit  which  was  proposed  to  solve  the  en- 
tire lateral  set  Is  shown  in  detail  In  Figure  II,   It  is,  however,  based 
on  the  same  ideas.  The  three  equations  which  are  satisfied  at  the  forward 
loop  amplifier  inputs  are: 

^recorded  "  ^g  "  ^gust  response  ~  "control  response 

^recorded  ~  ^gust  response  ~  ^control  response 

I'recorded  ~  H^gust  response  ~  ^control  response 

It  Is  to  be  noted  that  in  the  case  of  3,  the  recorded  quantity 

from  the  tape  was  based  on  an  accelerometer  measurement  and,  thus,  con- 
tained the  quantity  Pg,  In  addition  to  Pg^st  response'  "^^^   aircraft  re- 
sponse to  3g.  Therefore,  to  satisfy  the  above  equation,  the  Pg  was  fed 
back  to  the  amplifier  Input.   In  the  case  of  9  and  ift,  the  aircraft  response 
only  was  Included  in  tho  measured  \&?iponse   since  fhe   sensors  were   rate 
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gyros.  Approxrmatejy  Ijifinlte  resistance  feedback  was  therefore  provided 
around  the  forward  loop  amplifiers  in  these  two  channels. 

It  was  discovered,  however,  that  the  circuit  shown  in  Figure 
II  was  not  satisfactory.  During  operation  of  the  circurt  the  output  quan- 
tity -  Pg  exhibited  a  divergence^     This  Indicated  sonre  positive  feedback 
within  the  multiloop  system,   ft  can  be  seen  that  there  were  nine  loops 
through  the  aircraft  af^aLog^  ti   seemed   probable  that  the  -  3g  oufput   was 
controlled  primarily  by  the    loop  through  amplifier   9  which    is  shown  below 
schematical ly: 

— vvv^ 1  I WV^- 


^/Z 


I AA/VN- 


/5 


-'Wv'^ 


1 


This  circuit  obviously  provides  positive  feedback  and  is  thus 
unstable.  This  problem  emphasized  the  point  that  although  the  equations 
were  physically  satisfied  in  the  circuit  as  shown,  it  was  necessary  to 
Insure  system  stability  in  the  electronic  sense.  A  study  of  the  circuit 
showed  that  there  were  additional  unstable  loops,  but  that  their  effect 
was  minor  compared  to  the  effect  of  the  stable  loops  which  fed  back  to 
the  same  points. 

In  order  to  provide  stability,  the  analysis  circuit  was  modified 
to  the  form  shown  in  Figure  12.  The  difference  here  is  that  the  -  pg  quan- 
tity input  to  the  analog  is  now  taken  directly  from  the  recorded   data.  The 
analog  circuit  was  modified  to  maintain  satisfaction  oi  the   equations  of 
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motlon.  Circuit  stability  was  now  attained  since  the  unstable  -  3q  loop 
through  amplifier  9  was  now  eliminated. 

To  test  the  validity  of  the  analyzer  circuit,  a  preliminary  study 
was  made  in  which  the  aircraft  analogue  was  forced  by  a  ficticious  0g  quan- 
tity as  obtained  from  a  potentiometer  output.  The  aircraft  response  in 
roll  rate,  yaw  rate,  and  sideslip  was  recorded  on  tape.   It  is  to  be  noted 
that  the  recorded  sideslip  was  made  the  sum  of  the  analog  p  response  and 
the  forcing  function,  pg.  This  was  to  provide  similarity  to  the  actual 
flight  data  where,  as  previously  mentioned,  the  p  response  was  based  on 
an  accelerometer  measurement.  These  recorded  responses  were  then  played 
Into  the  analyzer  circuit  In  an  attempt  to  duplicate  the  original  forcing 
function. 

Figure  13  shows  the  original  forcing  function  and  the  recorded 
aircraft  responses.  When  these  responses  were  played  back  into  the  ana- 
lyzer circuit,  the  cpg  and  \g   outputs  were  essentially  zero,  as  expected. 
The  pg  output  was  found  to  be  a  good  reproduction  of  the  Input  wave  form, 
although  there  was  a  very  low  frequency  oscillation  superimposed  on  the 
trace.  This  is  shown  in  Figure  14(a),  where  the  solid  line  is  a  tracing 
of  the  original  p  input.   |t  was  determined  experimentally  that  this  low 
frequency  drift  was  strongly  affected  by  the  biasing  adjustment  of  the 
tape  output  from  the  ground  station  equipment.  Careful  adjustment  of  the 
bias  practically  eliminated  the  drift;  but  since  this  method  would  not  be 
feasible  with  actual  flight  test  data,  another  method  was  sought.   It  was 
decided  to  utilize  a  high  pass  filter  to  eliminate  the    low  frequency  os- 
cillation. This  was  connected  as  shown: 
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Since  the  use  of  such  a  filter  would  result  also  in  the  attenua- 
tion of  some  of  the  higher  frequencies,  the  corner  frequency  was  selected 
on  the  basis  of  compromise  between  minimum  attenuation  and  m\r\\m\jm   diver- 
gence. The  results  are  shown  in  Figure  14  (b),  (c),  and  (d).  A  corner 
frequency  of  1/44  cps,  as  provided  by  Rf  =  7  meg  Q,  seemed  to  give  the 
best  reproduction. 

Throughout  this  discussion  It  has  been  assumed  that  the  aircraft 
response  data  were  available  in  the  form  of  3,  cp,  and  ^,      It  can  be  seen 
from  Table  I  that  lateral  acceleration,  ay,  was  recorded  during  the  flight 
tests,  rather  than  sideslip.   It  was  necessary  to  obtain  3  from  this  data 
by  means  of  the  following  form  of  the  side  force  equation: 


(25)  +  3  = 


(p/2)Vo^S  C 


'^  TT  rJ~  ^y 


yp 


Cyp  2Vo    Cyp  2Vo    Cyp 


The  summing  and  scaling  circuit  utilized  to  solve  this  equation 
for  application  to  the  analyzer  is  shown  below; 
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Since  an  accurate  representation  of  p  depended  on  reliable  valuer 
°^  Cyo,  Cy  ,  Cy  ,  and  Cy.  ,  and  since  these  derivatives  were  not  Involved 
in  the  previous  matching  process,  it  was  necessary  to  check  their  validity. 
The  method  used  was  to  compare  the  quantity  ay  from  the  flight  data  with 
the  derived  ay  obtained  by  summing  the  appropriate  analog  outputs  in  re- 
sponse to  a  rudder  pulse.  The  results  are  shown  in  Figures  4,  5,  and  6, 
It  was  felt  that  the  matching  was  accurate  enough  to  justify  the  use  of 
P,  as  obtained  from  equation  25.  The  resulting  equation  for  ay  Is: 


(26)  +  a 


y  =  (P/2)Vn^S  Cw^  6+(p/2)Vo^S  b  C^  yH-  (P/2)Vo^S  b  Cw  i 
m     "^P       m    2Vo   P       m    Xio 


+  (P/2)^o  S  Cy^  6 
m       f 
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The  summing  and  scaling  circuit  for  obtaining  ay  for  the  matching 
process  is  shown  belov;: 
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The  analyzer  circuit  designed  for  use  with  the  JongitudinaJ  equa- 
tions was  based  on  the  same  principles  which  governed  the  lateral  circuit. 
The  schematic  of  this  circuit  is  shown  in  F'gure  15,  The  following  equa- 
tions were  satisfied  at  the  inputs  to  the  forward  loop  amplifiers: 

•  •  • 

^aircraft  response  ~  ®gust  response  ■*■  ^control  response 

Aircraft  response  ^  ^-gust  **■  ^-gust  response  "^  °^control  response 

It  should  be  noted  that  the  reliability  of  the  analyzer  output 
was  prpbabjy  improved  if  the  aircraft  maintained  a  very  steady,  smooth 
flight  path  during  the  carrier  approaches.  The  reason  for  this  is  that 
with  the  aircraft  responses  very  small,  the  aircraft  analogue,  in  the  ana- 
lyzer feedback  loops^  was  relatively  inactive.  This  was  henef/ciai  since 
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the  analogue  had  known  Inaccuracies.  It  was  therefore  desirable  that  the 
aircraft  stability  augmentation  equipment  be  operative  as  an  aid  In  main- 
taining steady  flight  during  the  carrier  approaches. 


-36- 


IV  RESULTS 

The  next  step  fn  the  data  analysis  was  to  play  Into  the  analyzer 
circuits  the  appropriate  quantities  from  the  flight  test  tape  recordings. 
Of  the  five  flight  test  tapes,  one  was  determined  to  be  virtually  unusable 
due  to  superimposed  noise.  The  four  remaining  tapes  were  each  played  into 
the  analyzer  and  Sanborn  recordings  were  made  of  the  five  gust  quantities 
for  each  carrier  approach.  Also  recorded  for  reference  purposes  were   the 
aircraft  response  in  five  degrees  of  freedom   and  the  various  control 
def lections. 

Figures  16  through  31  are  reproductions  of  gust  quantities  from 
the  lateral  and  longitudinal  analyzers.  Also  shown  are  traces  of  aircraft 
response  and  control  deflection.  Aircraft  altitude  also  appears  on  the  re- 
cordings as  an  aid  in  orientation.   It  was  apparent  from  the  traces  that 
some  drift  was  still  present  In  the  pg  channel,  probably  due  to  incorrect 
biasing  at  the  ground  station  output.  Further  Increase  of  the  filter  cor- 
ner frequency  was  not  practical  due  to  the  extensive  attenuation  of  the 
higher  frequency  components. 

Those  approaches  which  are  missing  from  the  set  of  Sanborn  re- 
cordings were  unusable  due  to  the  effects  of  "dropouts"  during  the  analyzer 
operation.  As  mentioned  previously,  the  tape  of  the  second  flight,  con- 
taining approaches  8  through  13,  was  completely  unusable  due  to  poor  quality 
of  recorded  information. 

Study  of  the  traces  revealed  a  number  of  points.   It  was  expec- 
ted that  if  the  analyzer  circuits  were  completely  valid  the  output  gust 
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quantltles  would  be  completely  random  In  character,  reflecting  the  nature 
of  the  actual  turbulence.  However,  Figure  19,  for  example,  shows  that 
there  were  occasional  sections  of  the  various  gust  traces  which  were  perio- 
dic in  character.  Conversely,  many  of  the  remaining  traces.  Figure  16  for 
example,  appear  quite  random.  Some  of  the  oscillations  appearing  were  of 
the  approximate  character  of  the  aircraft  short  period  and  Dutch  roll  oc- 
ci  Nation.  The  implication  here  was  that  the  FAB   analogue  was  not  com- 
pletely accurate.  As  a  result,  some  of  the  aircraft  response  characteris- 
tics appeared  in  the  analyzer  output. 

In  spite  of  the  apparent  imperfections  in  the  turbulence  data, 
some  qualitative  information  was  immediately  available.  The  records  showed 
that  for  all  approaches  there  was  a  marked  increase  in  the  magnitude  of  the 
gust  quantities  in  the  few  seconds  prior  to  aircraft  touchdown.  This  was 
a  verification  of  the  existence  of  the  turbulent  "burble"  immediately  astern 
of  an  aircraft  carrier  deck. 

In  order  to  detect  any  additional  trends  in  the  turbulence  quan- 
itites,  the  RMS  value  of  each  quantity  for  each  approach  was  determined 
with  the  aid  of  the  following  circuit,  which  included  the  Donner   Model 
3751  Function  Generator: 
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(27)     RMS  =    \  -1-  f  '<(ein^)cit 

I  Kt  J 
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In  order  to  eliminate  the  effect  of  any  DC  bias  or  low  frequency 
oscillations,  the  various  gust  quantities  were  channeled  through  a  high 
pass  filter  prior  to  entry  into  the  above  circuit.  A  corner  frequency  of 
,20  cps  was  selected  for  the  filter,  since  it  was  felt  that  frequencies 
below  this  would  probably  not  have  a  significant  effect  on  pilot  opinion. 
The  RMS  values  of  the  various  quantities  were  tabulated  and  are 
shown  in  Table  II.  Omissions  in  the  Table  were  due  to  the  effects  of  tele- 
metering  dropouts.  These  tabulated  values  do  not  include  the  effect  of  the 
final  "burble".  No  apparent  trends  were  determined  from  study  of  the  RMS 
tabulation,  except  for  an  indication  of  the  general  level  of  each  quantity. 
These  levels  remained  approximately  constant  during  the  course  of  the  ap- 
proaches, indicating  that  atmospheric  conditions  also  generally  remained 
the  same. 

It  was  also  noted,  with  reference  to  Table  III  that  the  angular 
velocity  gust  quantities,  cpg,  if^,    and  9^  were  all  of  approximately  the 
same  magnitude.  The  angular  quantities  3g  and  an  were  also  approximately 
equal.  This  was  a  strong  indication  that  the  turbulence  away  from  the 
Immediate  vicinity  of  the  ship  was  indeed  isotropic,  as  is  normally  assumed. 

It  was  expected  that  the  magnitude  of  the  gust  quantities  within 
the  "burble"  would  be  a  reflection  of  the  existing  relative  wind  conditions 
over  the  carrier  deck.  Although  no  formal  tabulation  was  made  of  the  RMS 
values  of  the  gust  quantities  within  the  "burble",  study  of  the  traces 
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revealed  that  these  RMS  values  were  approximately  doubled  In  the  final 
three  seconds  of  each  approach.  There  was  no  apparent  difference  in  the 
character  of  the  "burbles"  with  varying  relative  wind  conditions. 

It  was  noted,  however,  from  study  of  the  gust  traces  that  the 
"burble"  was  Initially  encountered  at  varying  distances  from  the  ship. 
Figure  32  shows  the  point  of  entering  the  "burble"  for  all  of  the  approaches 
represented  by  Figures  16  through  31,   It  is  interesting  to  note  that  those 
points  furthest  from  the  carrier  were  In  the  presence  of  high  relative 
wind,  with  no  x-wind  component.  Those  points  nearer  the  ship  were  with 
less  wind  velocity.  Finally,  those  points  appearing  2  and  3  seconds  prior 
to  touchdown  were  In  the  presence  of  "right  to  left"  x-wind,  indicating 
that  the  turbulence  trail  was  being  blown  across  the  flightpath  at  this 
point. 

In  Figure  33  are  plotted  the  pertinent  lateral  mode  data  from  a 
typical  approach.  The  increase  in  the  magnitude  of  the  gust  quantities  in 
the  final  few  seconds  before  touchdown  is  apparent.  However,  it  was  also 
noted  that  the  gust  quantities  at  this  point  had  some  periodic  characteris- 
tics.  In  addition,  the  oscillations  are  approximately  of  the  same  period 
as  the  aircraft  roll  and  yaw  response.   In  fact.  In  the  case  of  cp  and  9g, 
the  oscillations  are  exactly  180°  out  of  phase.  Because  of  apparent  "im- 
purities" such  as  this  in  the  turbulence  data  it  was  felt  that  a  more  de- 
tailed statistical  analysis  of  the  data  should  be  made  to  determine  its 
true  spectral  character.  Such  an   analysis  is  beyond   the  present  scope 
of  thi s  project. 
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V  CONCLUSIONS  AND  RECOMMENDATIONS 

Based  on  the  results  of  this  Investigation,  it  is  concluded  that: 
I.  Accurate  representations  of  the  atmospheric  turbulence  may  be  obtained 
by  the  analysis  method  described  in  the  report,  assuming  that  the  aircraft 
analogue  is  val id. 

2.  The  consistent  similarity  in  the  magnitude  of  the  angular  velocity 
gust  quantities,  and  the  similarity  in  the  magnitudes  of  gust  angular  quan- 
tities, as  shown  in  Table  IV,  indicated  that  the  turbulence  away  from  the 
Immediate  vicinity  of  the  carrier  was  isotropic  in  character. 

3,  The  atmospheric  turbulence  in  the  vicinity  of  an  aircraft  carrier  is 
characterized  by  a  strong  "burble"  In  the  vicinity  of  the  approach  end  of 
the  carrier  deck.   It  is  not  possible,  based  on  these  results,  to  give 

a  quantitative  description  of  this  turbulence  although  it  appeared  that  t/ie 
character  of  the  "burble"did  not  vary  noticeably  with  changing  relative 
wind  conditions.  However,  the  point  at  which  the  "burble"  is  initially 
encountered  is  a  function  of  the  relative  wind  conditions. 
It  is  recommended  that: 

1.  A  detai  led  spectrum  analysis  be  made  of  the  turbulence  quantities  re- 
sulting from  this  investigation,  in  order  that  undesirable  portions  of  the 
spectrum  be  eliminated,  making  the  information  suitable  for  application 

to  the  variable  stability  Navion. 

2.  A  further  statistical  analysis  of  this  data  be  made  to  determine  quan- 
titatively the  character  of  its  components. 

3.  An  investigation  be  conducted  into  the  applicability  of  this  analysis  cir- 
cuit for  the  general  solution  of  complex  Inverse  transfer  function  problems. 
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TABLE  I 
RECORDED  FLIGHT  TEST  DATA  QUANTITIES 


QUANTITY 
Rol I  rate,  e.g. 
Pitch  rate,  e.g. 
Yaw  rate,  e.g. 
Rol I  rate,  seat 

Pitch  rate,  seat 

Yaw  rate,  seat 


Pitch  and  rol I  rate, 
vertical  gyro 

Lateral  acceleration 


Normal  acceleration 
Rudder  deflection 
Left  ai  leron  def I . 
Left  spoi  ler  def I . 
Stabi  lator  def I  . 
Rudder  pedal  def I . 
Lateral  stick  def I , 


SYmOL 

«Pcg 

9cg 

• 

teg 

^'seat 

• 
®seat 

*seat 

• 

®VG 
ay 

az 
6r 

6,. 


Sped 
^stick,  I  at 


SENSOR 

Humphrey  rate  gyro.  Model  01-0238 

Humphrey  rate  gyro.  Model  01-0238 

Humphrey  rate  gyro.  Model  01-0238 

Minn.-Honeywei I  rate   gyro.  Model 
JG  7005A4-9 

Minn.-Honeywei J   rate  gyro.   Model 
JG  7005A4-9 

Minn. -Honeywell  rate  gyro.  Model 
JG  7005A4-9 

Minn. -Honeywell  rate  gyro.  Model 
JG  7044A4 

Conner  linear  accelerometer. 
Model  4310 

Donner  I i near  acce I erome ter. 
Model  4310 

2  K  0,  modified,  spring  loaded, 
selsyn  potentiometer 

IRC  2  K  Q  dual,  spring  loaded, 
potentiometer 

Bournes  4  inch  1 1  near 
potentiometer 

2  K  Q,  modified,  spring  loaded, 
selsyn  potentiometer 

spring  loaded,  modified  selsyn 
potentiometer 

spring  loaded,  modified  selsyn 
potentiometer 


i 


-A^,- 
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TABLE  I 
(Continued) 


QUANTITY 
Longitudinal  stick 


Left  and  right 
throttle  position 


Airspeed 
Altitude 


SYhCQL 


SENSOR 


^stick   long  spring  loaded,  modified   selsyn 
potentiometer 


'%'  "U 


h 


spring  loaded,  modified  selsyn 
potentiometers 

W.  0.  Leonard  Mbdel  50-2100 

W.  0.  Leonard  Wbdei  50-2000-9 
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TABLE  II 

RECORDED  FLIGHT  DATA  QUANTITIES, 

TELEMETERING  CHANNELS  UTILIZED 


QUANTITY 

Ro 1 1  rate,  eg 

Pitch  rate,  eg 

Yaw  rate,  eg 

Roll  rate,  seat. 

PIteh  rate,  seat. 

Yaw  rate,  seat, 

PIteh  rate,  vert  lea  I  gyro 

Roll  rate,  vertical  gyro 

Lateral  acceleration 

Normal  acceleration 

Rudder  deflection 

Left  aileron  deflection 

Left  spoiler  deflection 

Stab i  later  deflection 

Rudder  pedal  deflection 

Lateral  stick  deflection 

Longitudinal  stick  deflection 

Left  throttle 

Right  throttle 

Airspeed 

Altitude 


CHANNELS  UTILIZED 
3 
3 
3 
3 
3 
3 
I 
I 

2 
2 
2 
2 
2 
2 
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TABLE  III 
WIND  CONDITIONS  DURING  CARRIER  TEST  APPROACHES 
(Headwind  and  Crosswind  Components) 


RELATIVE 

RELATIVE 

APPROACH  mmFR 

HEADWIND  (KNOTS) 

CROSS];JlND   (KNOTSl 

1 

36 

0 

2 

36 

0 

3 

36 

0 

4 

36 

0 

5 

36 

0 

6 

36 

0 

7 

36 

0 

8 

19 

0 

9 

19 

0 

10 

19 

0 

II 

29 

0 

12 

26 

0 

13 

37 

0 

14 

31 

0 

15 

28 

0 

16 

26 

0 

17 

25 

0 

18 

25 

0 

19 

25 

0 

20 

23 

0 

21 

23 

0 

22 

23 

0 

23 

23 

0 

24 

26 

0 

25 

24 

7  1/2  right 

26 

28 

5  right 

27 

24 

5  right 

28 

27 

5  right 

29 

35 

0 

30 

34 

6  right 

31 

34 

5  right 

32 

34 

9  right 

33 

30 

5  right 

APPROACH 
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TABLE  IV 
RMS  VALUES  OF  TURBULENCE  QUANTITIES 


,    ^  . 

Qg 

3g 

9g 

!g 

(deqrees) 

(deq/sfic) 

(decrees) 

(dec/ sec) 

(dea/secl 

2.21 

3.06 

1.94 

2.88 

3.39 

2.21 

2.78 

1.69 

- 

- 

2.39 

2.89 

1.47 

2,77 

4.18 

2.21 

3.41 

- 

2.89 

3.86 

2.16 

2.72 

1.75 

3.19 

3.55 

2.16 

2.77 

1.66 

3.27 

3.75 

2.16 

3.17 

DATA 

1.79 
UNUSABLE 

3.27 

3.73 

I 

2 

3 

4 

5 

6 

7 

8 

9 
10 
II 
12 
13 

14  1.64        -         1.86       2.70       3.76 

15  -         3.60       2.03       2.71       3.66 

16  1.66        -  -         3.34 

17  1.63       2.69        -  -  - 

18  1.53       3.17       1.97        -         3.11 

19  1.43       2.99       1.97       3.03       3.66 

20  -  -         1.56       3.34       4.25 

21  -         2.97        -         2.93       3.34 

22  2.32       2.97       1.44       3.33       4.24 

23  -         3.17        -         3.22       3.66 

24  2.31       3.14       1.61       3.51       3.89 

25  2.41       2.89       1.57       2.95       3.48 
25           2.09       3.34        -         3.04       3.78 

27  -         3.16        1.77        -         3.21 

28  -  »  _         -         - 

29  2.51        -         1.67       2.84 

30  2.23       3.17       1.98       3.14       3.76 

31  2.09       3.27       1.91       2.86       3.89 

32  2.43       3.17       2.10       2.99       3.67 

33  2.20       3.27       2.10       3.04       3.83 
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